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ABSTRACT 

We propose a functional form for the IMF, the IMF, which is a natural heavy-tailed ap- 
proximation to the log-normal distribution. It is composed of a low-mass power law and a 
high mass power-law which are smoothly joined together. Three parameters are needed to 
achieve this. The standard IMFs of Kroupa (2001, 2002) and Chabrier (2003a) (single stars 
or systems) are essentially indistinguishable from this form. Compared to other 3-parameter 
functions of the IMF, the IMF has the advantage that the cumulative distribution fimction 
and many other characteristic quantities have a closed form, the mass generating function, for 
example, can be written down explicitly. 

Key words: stars: luminosity function, mass function, methods: statistical, methods: data 
analysis 



1 INTRODUCTION 

The initial mass function of stars (IMF), the spectrum of stellar 
masses at their birth, is of fundamental importance in many fields 
of Astronomy. Since the seminal work of Salpeter (1955), who in- 
vestigated the power-law part of the massive stars, a huge observa- 
tional and theoretical effort has been made to constrain this distri- 
bution. Towards the lesser masses the IMF deviates from a power 
law and follows more a lognormal shape (Miller & Scalo 1979). At 
present, the whole shape of the IMF is usually described by power- 
law segments (Kroupa 2001, 2002) or by a lognormal segment plus 
a power law segment (Chabrier 2003a,b, 2005). The aim of this 
paper is to provide an alternative, practical functional form for the 
IMF together with all its characteristic quantities (see Table 1 for 
the formulae and Figs. 3 and 6) ^ . More observational and theoret- 
ical aspects of the IMF can be found in recent reviews (e.g. Scalo 
1986; Chabrier 2003a; Zinnecker & Yorke 2007; Elmegreen 2009; 
Bastian et al. 2010; Kroupa et al. 2011). 

The IMF is usually believed to be a smooth function over the 
whole mass range, from brown dwarfs to O stars. However, Thies 
& Kroupa (2007) and Thies & Kroupa (2008) argued that a sudden 
change in binarity properties around the hydrogen burning limit in- 
troduces a discontinuity in the single star IMF as well. This discon- 
tinuity in the single-star IMF can still lead to a system IMF without 
discontinuities over the whole mass range (Thies & Kroupa 2007; 
Kroupa et al. 2011). In view of the simplicity aspect of our pro- 
posed IMF form we neglect any discontinuity. 

The proposed functional form, the L3 IMF, fulfils several de- 
mands on the form of the IMF: It describes the whole (system) 
mass range with a single function. This has been achieved by sev- 
eral other functional forms as well (Larson 1998; Chabrier 2001; 
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Paresce & de Marchi 2000; Parravano et al. 2011; Cartwright & 
Whitworth 2012). However, compared to these forms, the L3 IMF 
has the advantage that its cumulative distribution function is in- 
vertible, so that sampling from the L3 IMF is very easy. No spe- 
cial functions (e.g. the error function) are involved to normalise 
the La IMF as a probability. Beyond that the analytical form al- 
lows also for simple, closed forms of characteristic quantities, such 
as the peak or the "breaks", the masses from which on the power 
laws reigns. Furthermore, with three parameters, two controlling 
the power-law behaviour at low and high masses and one location 
parameter, the number of parameters is as small as possible. 

The motivation for the L3 IMF is of purely pragmatic nature, 
it is a functional form that describes the data in a very practical way. 
It would be pleasing if the 1/3 IMF could be more "theoretically" 
motivated. One could try to find a connection to some generalised 
log-logistic growth processes, in analogy to logistic growth, as the 
La is related to the log-logistic distribution. However, it remains 
questionable whether such a (non-stochastic) growth theory would 
be capturing the star formation process in its entirety (cf the dis- 
cussion about logistic growth in Feller 1968, p. 52). Where would 
be the place of, for example, feedback or stellar dynamics in shap- 
ing the IMF if growth alone gives all parameters of the IMF? Thus 
it seems futile to follow such thoughts and we do not attempt to find 
any reasons for our proposed functional form, other than its utmost 
simplicity and practicality. 

The organisation of this paper is the following: After some 
general definitions we discuss in Section 2 established functional 
forms and required parameters of the IMF. The L3 and IMFs are 
motivated and defined in Section 3 as heavy-tailed approximations 
and extensions to the log-normal distribution. This is followed by a 
detailed description of the L3 IMF and its characteristic quantities 
in Section 4, the S4 IMF is discussed in Appendix A. Section 5 
gives the "canonical" parameters for the La IMF, matching it to the 
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Kroupa (2001, 2002) and Chabrier (2003a,b) IMFs. Sec. 6 contains 
the conclusions of this article. 



2 PROPERTIES OF THE IMF 
2.1 Definitions 

We normalise the IMF as a probahility density function (pdf), 
the IMF tells us about the relative frequencies of stars of various 
masses in linear mass space. This allows us to use common statis- 
tical techniques, e.g. to estimate the parameters. For functions nor- 
malised as pdf we use the symbol p{m), for their integrals, the cu- 
mulative distribution function, the symbol P(jn). The cumulative 
distribution function is related to the observed number frequency, 
N{m), by P{m) = ^N{m), where ritot is the total number of 
observed stars. The standard normalisation condition for a proba- 
bility is 



•I mi 



p{m)dm, 



(1) 



where mi and m„ are the lower and upper mass limit, respectively. 

Historically there exist two alternative descriptions of the IMF, 
in linear or in logarithmic space, the small-a and the big-F notation. 
The use of the IMF as probability of m leads naturally to the linear 
(small-a) description, the IMF is fulfils 



p{m) = 



dP{m) 

dm 



1 diV(m) 



ntot dm 



(2) 



A power law IMF has then the exponent —a, p{m) oc m^". In 
the logarithmic description the IMF is normalised as probability of 
log m, not m. 



piog(logm) = 



dPiog(logm) 
dlogm 



1 dNiog{logm) 



ntot 



dlogm 



(3) 



piog(log m) is connected to the linear pdf via 

, , dflog(logm) dlogm 1 , , 

dlogm dm m 

Thus, a power law pdf in m, p{m) oc m~", transforms into 
j3iog(m) cx m~^°'~^^ orpiog(m) oc m~^ , where F = a — 1. 

We define the exponent (sometimes referred to as "slope", but 
that should be reserved for the logarithmic description), as a func- 
tion of mass via 



S{m) 



dlogp(m) _ dlogp(r 
d log m dm 
A power-law IMF can then be written as 



p{m) 



-Sim) 



(5) 



(6) 



We follow the convention that the negative sign is not included in 
the exponent. Thus, in our notation the Salpeter (1955) exponent is 
positive, a =+2.35. 

2.2 Tlie standard IMFs and otiier functional forms 

The Kroupa (2001, 2002) single-star IMF consists only of power- 
law segments. 



PKroupa(m) = 



Akom-"-" 
Afeim-i-3 
Ak2m'^-^ 
Akam-^-^ 



0.01 M0< m< 0.08 Mq 
0.08 M0<m< 0.5 Mq 
0.5 Mq < m< 1 M© 
IM0 <m (< 150 M0) 



,(7) 



fcimj ^'^^^-^ and 



with feo = 1, ki 

kz = femg ^■^+^-^(= k-z) where mi = 0.08 Mq, m2 = 0.5 Mq 
and ma = 1 Mq (a practical algorithm for the calculation of the 
ki is given by Pflamm-Altenburg & Kroupa 2006). A is some 
global normalisation constant. This form is highly adaptable, which 
comes at the price of a large number of parameters. On the prac- 
tical side, the Kroupa (2001, 2002) IMF has tiie advantage that 
many derived quantities can be calculated without involving spe- 
cial functions (ctmiulative distribution function, quantile function, 
mean mass etc.), but with several "if" statements to specify the 
mass ranges. 

Chabrier (2003a,b) combined for the single-star IMF a log- 
normal distribution at the low-mass end with a high-mass power 

law. 



PChabrier(m) 



Aki^e H 
Ak2m-^-^ 



loglO 0.079 \2 



m < 1 Mq 
m > 1 Mq 



with fei = 0.158 and k2 = 0.0443 and the global normalisation 
constant A. The lognormal and the power-law part connect up more 
or less smoothly, without the "kinks" of several power-law seg- 
ments (although there is still the small kink at 1 Mq). Calculating 
the cumulative distribution function involves the error function, but 
random variates can be created without any specialised algorithms 
from standard Gaussian distributed random numbers. 

A piece-wise functional form of the IMF is somewhat unsatis- 
fying, and several alternatives covering the whole mass range have 
been proposed in the literature. There are, for example, the func- 
tional forms of Larson (1998) 

-("+!) 

PLarsona(m) OC ( 1 + — ) , (9) 



1 1+"^ 

m \ p 



and 



PLarsonb(m) OC m "e 

form 3 of Chabrier (2001), 

PChabrier3(m) OC m~"e 



(10) 



(11) 



or the tapered power law form of Paresce & de Marchi (2000), De 
Marchi et al. (2010), Hollenbach et al. (2005) and Parravano et al. 
(2011), 



PTapereredPL(m) OC m " ( 1 - e 



(12) 



The IMF forms of eqq. 9, 10, 1 1 and 12 are very similar to our pro- 
posed form of the IMF, but their integrals contain the incomplete 
gamma function or the hypergeometric function. A cumulative dis- 
tribution function without closed form is hard to invert, so that spe- 
cial algorithms are necessary for random variates from these distri- 
butions. 

Recently, Cartwright & Whitworth (2012) proposed a com- 
pletely different class of distribution functions for the IMF descrip- 
tion, stable distributions. Stable distribution (e.g., the Gaussian dis- 
tribution) arise naturally in the context of stochastic processes, of 
which the star formation process is one example. Related to stable 
distributions, and also the outcome of stochastic processes is the 
class of infinitely divisible distributions, such as the lognormal dis- 
tribution (e.g. Zinnecker 1984, Elmegreen & Mathieu 1983; Thorin 
1977 for infinite divisibility). The choice of stable distributions is 
motivated by their relation to stochastic processes, however, they 
are also used only as a fitting function, as the exact stochastic pro- 
cess describing star formation has not yet been formaUsed. Also, 
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Figure 1. Comparison of the Normal (dashed hne) and Logistic (solid line) 
probability density, with a logarithmic y axis. The Logistic distribution has 
heavier tails. 




X 



Figure 2. Comparison of the log-Normal (dashed Une) and log-Logistic 
(solid line) distribution, scaled to 1 at i; = 1. The dotted lines are at e~^, 
e~^,e^ and e^. The tails of the log-Logistic distribution are asymmetric. 



typically they do not have a closed form for the distribution func- 
tion itself, which is an important practical aspect. 



2.3 How many parameters for the IMF? 

The IMF seems to have a lognormal body with a power law tail 
on both the high-mass and the low-mass side. In order to describe 
this behaviour, four parameters appear to be required: a location 
parameter (which is not necessarily the "peak" or the mean), a 
scale or width parameter (which is not necessarily the variance), 
the low-mass and high-mass power-law exponents. There are no 
stars of zero or infinite mass, so that additionally an upper and a 
lower mass limit has to be introduced, so the total number of pa- 
rameters is 4-1-2. This is two parameters less than in the schematic 
IMF of Bastian et al. (2010), where additionally two "mass breaks" 
are introduced, i.e. 6-1-2 parameters. However, if one requires that 
the lognormal part merges smoothly into the power law tails, then 
the scale parameter sets the width of the IMF and consequently the 
mass breaks. The mass "breaks" are then not parameters any more, 
but derived quantities. 4-1-2 seem therefore to be the necessary nvrni- 
ber of parameters to describe the IMF. The B4, IMF discussed later 
is a smooth function over all masses and has the mentioned 4-1-2 
parameters. 

The number of parameters of the IMF can be reduced by one, 
because it is not necessary to explicitly include a scale parame- 
ter to fit the "canonical" IMF. Only a location parameter and the 
two exponents suffice to achieve this. Several 3-1-2 IMFs have been 
suggested in the literature (eq. 11, IMF 3 of Chabrier 2001; eq. 12 
Paresce & de Marchi 2000; De Marchi et al. 2010; HoUenbach et al. 
2005; Parravano et al. 2011). Our proposed L3 IMF also has only 
3-1-2 parameters. 

2-1-2 parameter functional forms (eq. 9 and 10) have been given 
by Larson (1998), with a location parameter and only a high-mass 
exponent. With only 2-1-2 parameters it is difficult to fit the low- 
mass end of the IMF. 

For comparison, the Kroupa (2001, 2002) has 5-1-2 parame- 
ters (three exponents, two thresholds, two limits) and the Chabrier 
(2003 a,b) IMF has 4-1-2 parameters (mean, variance, one exponent, 
one threshold, two limits). 



3 HEAVY-TAILED APPROXIMATIONS TO THE 
LOGNORMAL DISTRIBUTION 

Starting point for the search of a functional form for the IMF is the 
relation between the Normal distribution and the Logistic distribu- 
tion (see e.g. lohnson et al. 1994, 1995). The Normal distribution. 



1 



_1 (x-rtl 



2Tia 



(13) 



can be approximated in the central region for a = 1 by the Logistic 
distribution. 



PLix) 



(14) 



where cr' = e~ 2 . The ratio of the two probability densities is close 
to unity between ~2a and +2a, but drops off strongly outside. 
This behaviour is evident in a logarithmic plot of both densities 
(Fig. 1), the tails of the logistic distribution are much heavier than 
the normal distribution, with fixed exponents. 

In order to translate the relation of Normal and Logistic distri- 
bution to the lognormal distribution, 

1 



we rewrite the lognormal density function as 



1 -iltn 

PiogMlx) oc -e 

X 



(15) 



(16) 



Inserting In 
bution function, 

PlogL{x) = 



l + e" 



for into the Logistic cumulative distri- 



(17) 



and taking the derivative gives the log-Logistic density. 



piogz,(a;) oc 



(18) 



Figure 2 shows piogAT and piogi,, again with a = 1 and a' = 
e~ 2 . The log-logistic distribution follows the lognormal distribu- 
tion over about two orders of magnitude and deviates with asym- 
metric tails. 
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Table 1. Collection of formulae for the L3 form of the IMF. The values given for the parameters and characteristic quantities are to match the "canonical" 
single-star IMF (Kroupa 2001, 2002; Chabrier 2003a), values in parentheses for the "canonical" system (binary star) IMF (Chabrier 2003a). B{t; p, q) is the 
incomplete Beta function. For the limits we adopt the fiducial values mi = 0.01 M© and ruu = 150 Mq, which are only needed for the normalisation. 



(1) AuxiUiary function: 



G(m) : 



(2) 
(3) 



(5) 

(6) 
(7) 

(8) 



(9) 

(10) 
(11) 



Quantity 
Functional form 



Formula 



Cumulative distribution function (CDF) Pj^3 (m) 



Probability density function (pdf) 



(4) Quantile function 



Parameters 

High-mass exponent 

Low-mass exponent 

Scale parameter 

Lower mass limit 

Upper mass limit 

Shape characterising quantities 

Effective high-mass exponent 

Effective low-mass exponent 



Lower power-law mass limit 
Upper power-law mass limit 

Exponent (N.B. 5(oo) = a (= 2.3)) 

Scale characterising quantities 

Mean mass (Expectation value) 

Median mass 

Mode (most probable mass) 
"Peak" (maximum in log-log) 



G(m) - G{mi) 



G(m„) - G{mi) 

P„w^.(^)-"(..(^)' 



H G{mu) — G{mi) 



m{u) = fj, 



u(G{mu) - G{mifj + G{m,} 



^ 1 

1-/3 



a = 2.3 (2.3) 
13 = 1.4 (2.0) 
^J. = 0.2 (0.2) M0 
mi = 0.01 M0 
rrin = 150 M0 

a = 2.3 (2.3) 

-y = a + 13(1 - a) = 0.48 (-0.3) 

m-y = iJ,eT^ = 0.043 (0.043) Mq 
ma = lJ.e^^ = 0.93 (0.93) Mq 



S{m) = a + ^{1 - a) I 1+ ( — 



E{m) (= m) = expressed by Beta function, see eq. 25 
= 0.36 (0.62) Me 



m = fi 



(G{mu) - G(mi)^ + G(mi) 



0.10 (0.21) Mq 



I mi 7 > 

mp = //(/9- 1)S=T = 0.10 (0.20) Me 



Definition 



P{mi) = and P(mu) = 1 



dm 



m{u) = P-I(u),w e [0,1] 



a 1 (typically a > 0) 
13^1 (typically /3 > 0) 

> 
nil > 
mu > 

lim p(m) oc m~" 

lim p(m) oc m '' 
m— >0 

p{m 6 [mi . . . mj]) Ri m 
p{m e [ma ■ ■ ■ mu]) ~ m~" 
d log p(m) 



5(m) = 



d log TO 



E{m) = I mp{m)dn 
J mi 



fh = arg maxp(m) 

m 

dlog {mp{m)) 



d log TO 







The log-Logistic distribution of Fig. 2 already looks very 
much like the IMF. Only the high-mass and low-mass exponents are 
still fixed. In fact, this is not quite correct, because the meaning of 
cr' has been changed from the width of the distribution (i.e. a scale 
parameter) to determining the low-mass exponent (i.e. a shape pa- 
rameter). Arbitrary exponents for the low-mass and the high-mass 
tail can be introduced by writing 



PLsim) oc 




Unfortunately, j3 is not the exponent at low masses, which is the 
price paid for eq. 19 having a very simple cumulative distribution 
function. Probability densities similar to eq. 19 (two exponents and 



fj.) are known under several other names, particularly in economics. 
We will refer to it as generalised log-Logistic distribution, or in 
short "1/3 IMF", because it has three (shape) parameters. 

A parameter that changes the width of the IMF can be intro- 
duced by writing 



PB4(m) = ''^^t^+g) - (20) 

a is now the scale parameter, and q and ,5 the exponents of the 
power-law tails. The integral of eq. 20 does not have a closed form, 
but can be transformed to the incomplete Beta function. Therefore, 
probability densities of the type of eq. 20 are known as (gener- 
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Figure 3. Probability density function for the Lz functional form of the IMF 
(Table 1 eq. 3) with the "canonical" parameters, a = 2.3, /3 = 1.4, = 
0.2 and the limits m; = 0.01 M0 and rriu = 150 Mq. It follows approx- 
imately m-" for m > m„(Ri 0.9 Mq) and m-~< = m-("+''(l-")) 
for m < m-|,(Ri 0.04 Mq) with 7 = 0.48. Also shown are the locations 
of mean, median and mode, which are all different because of the skewed 
distribution. The infamous "peak" (maximum in log-log) is not the location 
at which the two power laws cross over. This happens at the scale parameter 



alised) Beta distributions. Because of the four parameters we will 
refer to it as B4, IMF. 

The following Sections will show, that the "canonical" IMF 
(Kroupa 2001, 2002; Chabrier 2003a) can be very satisfyingly de- 
scribed by the Li IMF. The introduction of a as an additional scale 
parameter seems not to be necessary. Therefore we consider in the 
following only the L3 IMF and give the corresponding equations 
and parameter values for the B4 IMF in appendix A. 



4 THE Lz IMF 
4.1 Functional form 

The probability density of the Lz IMF is given in eq. 19, or, with 
the normalisation constant, in Table 1 eq. 3. Table 1 collects all for- 
mulae for the Lz IMF. Figure 3 shows the Lz IMF with its charac- 
teristic quantities for the "canonical parameters" of the single-star 
IMF. The particular advantage of the Lz IMF is that the integral of 
the probability density is very simple. 



1+ 



dm oc 1 + 



= : G(m). (21) 



The full cumulative distribution function, including the upper and 
lower limits {mi and m„), is then 



P(m) 



G(m) - G{mi) 
G{mu) - G{mi) 



(22) 



(also eq. 2, Table 1). Eq. 22 can be readily inverted to give the quan- 
tile function (Eq. 4, Table 1). Generating a random mass from the 
Lz IMF (i.e. inserting a uniform random number u in the quantile 
function) can then essentially be done in a single line of code. 

The two shape parameters have different meanings for the L3 
IMF. For large masses limm^ 00 pim) oc m~°, i.e. a is the high- 
mass exponent. In order that the Lz IMF is defined a 7^ 1 is re- 




Figure 4. Q/37 plot, showing the value of the low-mass exponent (p(m) <x 
m~^. Table 1 eq. 5) as a function of a and /3. The lines solid for integer a 
and a + 1/2 and dashed for a + 1/4 and a + 3/4. The red dotted line is 
for a = 2.35. 
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Figure 5. ("Alpha plot") Exponent of the L3 IMF (black solid curve) and 
its approximation in log space (red solid lines). The points ma and m-y are 
defined as the intersection of the straight line approximation of the exponent 
at with the limiting exponents 7 and a. For comparison the Kroupa (2001, 
2002) IMF (blue dashed line) and the Chabrier (2003a) IMF (green dotted 
line) are given as well. 



quired, typically will be a > 1. For small masses the limiting case 
is lim„i^o p(?7i) oc with 7 = a + — a). Therefore the 
parameter /? is not the low-mass exponent. This inconvenience of 
/? and 7 is the trade-off for the very simple cumulative distribution. 
Again, in order for the Lz IMF to be defined /3 7^ 1 is required, 
typically will be /3 > 1. For a > 1 and /3 > 1 the largest value 
that 7 can take is +1, i.e. p(m) oc mT^ . 7 will be negative for 
/? > -^^zri - A graphical representation of the relation between the 
exponents is given in the "q/37 plot". Fig. 4, where the value of 7 
for given a and /? can easily be read off. 



4.2 Breakpoints 

Related to the low- and high-mass exponents is the question of the 
"breakpoints" in the IMF. As for the Lz (and the B^) IMF there 
is a smooth transition between the exponents, proper breakpoints 
do not exist. Nevertheless, it is useful to know from where the Lz 
IMF can be approximated by a power law. Our approach to find 
the breakpoints is via the exponent as function of mass (defined in 
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eq. 5, given for the Lz IMF in eq. 8, Table 1) For the L3 IMF the 
curve of the exponent vs. logm is "S"-shaped, see the black solid 
line in Fig. 5. This "S" shape can be approximated by three straight 
lines (red in Fig. 5), of which two are horizontal at 7 and a. The 
intermediate, increasing part follows 



AS{r, 



dlog m 



log 











(23) 



a straight line in logm. We define now the breakpoints, and 
ma, as the points where g{m^) = 7 and g{mc,) = a. Formulae 
are given in Table 1, eqq. 6 and 7. The agreement of the L3 IMF 
and the power-law segments below and above rria is good, as 
can be seen in Fig. 3, where the power-law segments are shown as 
red lines, which are in fact barely visible. 



4.3 Characteristic masses 

Characteristic mass scales of the L3 IMF are also shown in Fig. 3 
and given in Table I. Because the IMF is skewed, the mean, median 
(Eq. 9, Table 1) and mode (most probable value, eq. 10, Table 1) 
are all different. Also, note that ^ is not directly related to any of 
them, it is the inflexion point of the exponent. Calculating the mean 
of the L3 IMF involves incomplete Beta functions^ (B{x;p, q) — 
jx jp-ij-j^ _ ty^^dt). Using the transformation 



1+ 



t{m) 



the mean can be expressed as 

B (t(m„); a, b) 



E{m) 



B {t{mi); a, b) 



G{mu) — G{rni) 



(24) 



(25) 



where a = and b — j3 — and G{m) is the auxiliary 
function given in Table 1 , eq. 1 . 

The "peak" of the IMF refers to the maximum in the logarith- 
mic description. The also very simple formula for mp is given in 
eq. 11, Table 1). 



5 "CANONICAL" PARAMETERS FOR THE L3 IMF 

Observationally, the shape of the IMF is constrained mainly by the 
number ratios of different mass ranges to each other, for example 
the ratio of high-mass to low-mass stars. Thus, a first approach to 
find the "canonical" parameters for the L3 IMF could be a fit to 
the cumulative distributions of the Kroupa or Chabrier IMF. This 
could be done in some objective way, for example by matching 
histograms of L3 to Kroupa or Chabrier. However, there are more 
properties that a "canonically" parametrised IMF should fulfil: Not 
only the number ratios, but also the mass ratios, the shape and the 
exponent should agree with each other We could not find an "ob- 
jective" procedure that would fit these constraints such that for all 



The incomplete Beta function is available in many scripting languages 
for data processing (R (open source), IDL etc.) or via Numerical recipes 
(Press et al. 2007). Sometimes what is called "incomplete Beta func- 
tion" is actually the regularised incomplete Beta function, Ix {p, q) = 
B{x;p, q) / B{p,q). This is the case for the functions PBETA in R and 
IBETA in IDL. 
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Figure 6. Probability density function for the L3 IMF, like Figure 3 but 
using the logarithmic description, shown together with the Kroupa (2001, 
2002) IMF (dashed blue line) and the Chabrier (2003a) IMF (green dotted 
line). 



of them the fit is good, the high-mass power-law tail leads to prob- 
lems. Therefore we choose the parameters "by hand" for an optimal 
agreement of the Lz with Kroupa and Chabrier in all the criteria. 

For observational data objective fits are, of course, possible, 
for example with the maximum likelihood method. There not only 
the upper mass exponent and the lower-mass exponent, but also the 
scale parameter p can be estimated. This is an advantage compared 
to the piecewise defined IMFs, where typically the "breakpoints" 
are not estimated. It is also possible to estimate the limits, in partic- 
ular rriu, which can also vary between star forming regions (cf. e.g. 
Weidner & Kroupa 2006, Maschberger & Clarke 2008, or Weidner 
et al. 2010 for an observational perspective and Maschberger et al. 
2010 for a varying rUu in simulations). 

In order to normalise the IMFs to be able to find the "canoni- 
cal" we choose mi = 0.01 Mq, near the deuterium burning limit. 
We set rrii, = 150 Mq, as this is commonly assumed (cf. Weidner 
& Kroupa 2004; Oey & Clarke 2005; Figer 2005), but are aware 
that in some star forming regions m^ can be at much higher masses 
(Crowther et al. 2010). As m„ lies well in the power-law tail, the 
exact value of it does not affect the parameter determination, a, /3 
and fj, are mainly constrained by the behaviour of the IMF below 
ma. 



5.1 Lz single star IMF 

In Fig. 6 we show in the logarithmic description the Lz IMF with 
parameters chosen such that it fits the "canonical" single-star IMF 
(a = 2.3, P — 1.4 and /i = 0.2 M0). For comparison we also 
show the Kroupa (2001, 2002) IMF and the Chabrier (2003a) IMF, 
both also normalised as probabilites. The difference between L3 
and Chabrier (2003a) is marginal, between Lz and Kroupa (2001, 
2002) equal to the difference between Kroupa (2001, 2002) and 
Chabrier (2003a). The effective low-mass exponent is 7 = 0.48 
for m < m-y = 0.042 M©. The high-mass break occurs at = 
0.93 M0, comparable to the start of the high mass power law of 
Chabrier (2003a) at 1 M©. In the Kroupa (2001, 2002) IMF the 
high-mass power law continues to 0.5 Mq . 

A comparison of the number fraction of stars in several mass 
bins is shown in the top panel of Fig. 7. The agreement between Lz 
and Chabrier (2003a) is again very good. The fraction of stars in 
the mass range 0.6-2 M© is slightly smaller for Lz, because of the 
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Figure 7. Comparison of the fractions of number (upper panel) and frac- 
tions of mass (lower panel) for the L3 IMF (solid), Kroupa (2001, 2002) 
IMF (blue dashed) and Chabrier (2003a) IMF (green dotted). The L3 IMF 
agrees better with the Chabrier (2003a) IMF, except for the mass range 
around I M0, where the power law is mounted onto the lognormal in the 
Chabrier (2003a) IMF 



smooth transition to the high-mass power law. There are differences 
between L3 and the Kroupa (2001, 2002) IMF at O.S-IM© and at 
0.01 - 0.3 M0, caused by the segments in the Kroupa form. The 
lower panel of Fig. 7 shows the fraction of total mass in the mass 
bins, 



%m 



100 



J^'' mp{m)&m 



(26) 



{nia and mt being the bin limits). La again agrees very well with 
Chabrier (2003a) and well with Kroupa (2001, 2002). 

As a last point we compare the exponent of the L3 IMF with 
Kroupa (2001, 2002) and Chabrier (2003a), see Fig. 5. Interest- 
ingly, although the probability density function, the cumulative dis- 
tribution function (fraction of stars, top panel of Fig. 7) and the 
mass distribution function (fraction of mass, bottom panel of Fig. 
7) of the L3 IMF agree more with a Chabrier (2003a), the exponent 
of the L3 IMF follows more closely the Kroupa (2001, 2002) IMF 



1 


1 


1 — 1 1 — 1 1 — r-TT 




1 


0.1 


_mi 


^ ; \ m mp ' 1 

^ m X. 

m nij nia % 




: 


0.01 










10-3 










10-" 


- 1 










0.01 


0.1 1 

nlSystem [Mg,] 


10 


100 



Figure 8. Comparison of the L3 system mass function (solid) and the 
Chabrier (2003a) system mass function (green dotted) in the logarithmic 
description. 
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Figure 9. Comparison of the fractions of number (upper panel) and frac- 
tions of mass (lower panel) for the L3 system IMF (solid) and Chabiier 
(2003a) system IMF (green dotted). As for the single-star IMF (Fig. 7) the 
lognormal-power law transition of the Chabrier (2003a) IMF around 1 Mq 
cannot be fitted exactly. 



5.2 L3 system IMF 

The system IMF for m < 1 has been given by Chabrier 
(2003a), 

PCtabner2003.System m = ylO.086 — e H 0.57 ) ^ (27) 

m 

and, with slightly modified parameters by Chabrier (2005), 

1 1 / logior,.-logio0.25 N2 

PCtabne,-2005.Sys.em(m) = A0.076 — e ^\ ) . (28) 

m 

{A is a normalisation constant). Above 1 Mq the system IMF fol- 
lows a power law with exponent 2.35 both in Chabrier (2003a) 
and Chabrier (2005). We adopt the Chabrier (2003a) form for 
m < 1 Mq and a power law with exponent 2.3 

The best parameters for Lg to fit the Chabrier (2003a) system 
IMF are a = 2.3, /3 = 2 and = 0.2 M©, taking m; = 0.01 Mq 
and m„ = 150 M©. A graph of both IMFs in the logarithmic 
description is given in Figure 8, where very good agreement is 
achieved. 



The effective low-mass exponent is then 7 = —0.3 with 
breakpoint m-y — 0.043 Mq and high-mass breakpoint at nia = 
0.93 Mq. The mean mass is 0.62 Mq which compares well with 
the 0.64 Mq for the Chabrier (2003a) system IMF. The mass for 
the median (0.21 Mq), the "peak" (0.20 Mq) and the mass scale 
parameter (/i — 0.20 Mq), by chance, coincide. Another coinci- 
dence is the near-equality of m-y and the mode (m — 0.042 Mq). 

As for the single star IMF, the fraction of stars and the frac- 
tions of mass over the range of mass bins is very comparable for 
the 1/3 system IMF and the Chabrier (2003a) system IMF (Figure 
9). 



6 SUMMARY 

The Lg IMF, a functional form of the IMF generalising the log- 
Logistic distribution, describes the whole stellar mass range with a 
minimum number of parameters (3 shape, 2 limits, see Table 1 that 
collects all formulae). It consists of a low-mass and a high-mass 
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power law that are joined smoothly together. Due to its analyti- 
cal simplicity many characteristic quantities (e.g. peak and mass 
breaks) can be given explicitly. The cumulative distribution func- 
tion is analytically invertible, so that drawing random masses from 
the L'i IMF is also very simple and does not involve a large pro- 
gramming effort. 

We have determined the parameters that fit the Lz IMF to the 
widely used single-star IMFs of Kroupa (2001, 2002) and Chabrier 
(2003a) and the system IMF of Chabrier (2003a). The L3 IMF fol- 
lows these IMFs very well, obtaining the same number and mass 
fractions of various mass ranges, so that it is an viable alternative 
functional form. 
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APPENDIX A: THE B4 IMF 

In some cases it might be necessary to include explicitly the width 
of the IMF in its functional form. For this an additional param- 
eter has to be introduced, so that the total nvunber of parameters 
is 4-1-2, two exponents, one location and one scale parameter, plus 
the obligatory upper and lower mass limit. The Li IMF can be 
extended to include this additional parameter, the cumulative dis- 
tribution function then contains then Beta functions, thus the name 
Bi IMF. We give the formulae for the B4. IMF in Table Al. The 
"canonical" parameters were determined "manually", as for the L3 
IMF, and the agreement with the IMFs for single stars (Kroupa 
2001, 2002; Chabrier 2003a) and the IMF for systems (Chabrier 
2003a) is comparably good. Figures Al and A2 show this in the 
logarithmic description. 

This paper has been typeset from a TgXy file prepared by the 
author. 
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Table Al. Collection of formulae for the B4 form of the IMF. The values given for the pai'ameter are to match the "canonical" single-star IMF (Kroupa 2001, 
2002; Chabrier 2003a), values in parentheses for the "canonical" system (binary star) IMF (Chabrier 2003a). B(t; p, q) is the incomplete Beta function. 



(1) Auxilliary function: 



Bi{m) = B {t{m); (t(/3 + i),(7{a - i)) t(m) 



Ml 



Quantity 
Functional form 



Formula 



Definition 



(2) Cumulative distribution function (CDF) Pb4(»ti) 



Bl(m)-Bi(m;) 



Bi{mu) - Bi{mi) 



(3) Probability density function (pdf) ps4(m) = A- 



CT(a + l3) 



1 



1 



o-M Bi{m„) - Bi{mi) 



Parameters 

High-mass exponent 
Low-mass exponent 
Location Parameter 
Scale parameter 
Lower mass limit 
Upper mass limit 

Shape characterising quantities 

(4) Lower power-law mass limit 

(5) Upper power-law mass limit 

(6) Slope (N.B. 5(00) = +0 (= 2.35)) 

Scale characterising quantities 

(7) Mean mass (Expectation value) 

(8) Median mass 

(9) Mode (most probable mass) 

(10) "Peak" (maximum in log-log) mp = fj, 



a = 2.3 (2.3) 
^ = -0.15(0.4) 
^ = 0.15(0.20) M0 
o- = 0.85 (0.80) 
mi = 0.01 M0 
rrtu = 150 M0 

met = fie^" 

S{m) = -/3 + (a + /3)- 



1 + 



T^t X/ — ^ B2{m^) - B2{mi) 

E(m) (= mj = fi^^ 

Bi{mu) ~ Bi{mi) 

m = P~^(i) (no closed form) 
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Figure Al. Logaiithmic description of the B4 form of the IMF, like Fig. 
6 for L3. The blue dashed curve is the Ki'oupa (2001, 2002) IMF and the 
green dotted curve is the Chabrier (2003a) IMF for comparison. 
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Figure A2. Comparison of the B4 system mass function (solid) and the 
Chabrier (2003a) system mass function (green dashed) in the logarithmic 
description, like Fig. 8. 



